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Recently, it has been shown, that the pair density of the homogeneous electron gas 
can be parametrized in terms of 2-body wave functions (geminals), which are scat- 
tering solutions of an effective 2-body Schrodinger equation. For the corresponding 
scattering phase shifts, new sum rules are reported in this paper. These sum rules 
describe not only the normalization of the pair density (similar to the Friedel sum 
rule of solid state theory), but also the contraction of the 2-body reduced density 
matrix. This allows one to calculate also the momentum distribution, provided that 
the geminals are known from an appropriate screening of the Coulomb repulsion. 
An analysis is presented leading from the definitions and (contraction and spectral) 
properties of reduced density matrices to the Kimball- Overhauser approach and its 
generalizations. Thereby cumulants are used. Their size-extensivity is related to the 
thermodynamic limit. 
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Introduction 



The description of an electronic many-body system in its ground-state (GS) by means 
of 1-body orbitals if> K (l) with the short-hand 1 = (rj, a\) is well-known and widely used. 
For example in the configuration- interaction (CI) method \l, 3, ^] these orbitals are used 
to build up the iV-electron Slater determinants. The linearly combined Slater determinants 
form the GS-wave function . . . ,N) such that the electron density is given by p(l) = 
S K ^kIV'kWI 2 with a non-idempotent occupancy v K being between and 1 as a consequence 
of the Pauli principle, ^ K v K = N. In the Hartree-Fock approximation and in the Kohn- 
Sham treatment of the density-functional theory (DFT), cf. e.g. j^, the orbitals t/> k (1) 
are solutions of effective 1-body Schrodinger equations and the occupanies are idempotent, 
vt — Vk., using thereby the aufbau principle. The DFT shows the 'power and charme' of the 1- 
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body density p(l) p. In the theory of reduced density matrices (RDMs) |fil I7L l8l. I9L 11(1 ll lL Il2| 
and in the density matrix functional theory (DMFT) j], 3] the ip K (l) diagonalize the 1- 
body reduced density matrix (1-matrix) so that 7(1|1') = ^« ^(iV/tV^l') holds, what is 
called spectral resolution. The v K are again non-idempotent as a consequence of the electron 
correlation phenomenon such that the non-idempotency of 7 or v K can be used to measure 
*. « - — . ,1 < ^ < ^ flfldQ. Ko, 80MS th e 

u K describe the occupation band structure. The ^ K (l) are the natural orbitals and the k are 

1- body quantum numbers. 

The 1-matrix 7(1 11') is a 2-point function, not to be confused with another important 

2- point function, the pair density (PD) p2(l, 2). This PD shows its power and charm in the 
Fermi hole for parallel-spin pairs (due to the Pauli 'repulsion' and modified by the Coulomb 
repulsion), in the Coulomb hole for antiparallel-spin pairs (due to the Coulomb repulsion), in 
the coalescing (or on-top) cusp and curvature properties, and in the possibility to calculate 
and discuss particle number fluctuations in partial regions of the system (e.g. Daudel loges or 
"stockholder" loges or Bader basins or Wigner-Seitz cells or ■ • • ) as another (more sensitive) 
measure of the correlation strength with the conclusion '(strong) correlations (strongly) 
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suppress such fluctuations', cf. |1|], p. 157, and |17|, [18||- Strong electron correlation thus 



localizes electrons. Examples are the Wigner crystallization of the low-density HEG, the 
Wigner-like charge ordering in Yb4As3 19j |. and the insulating GS of CoO (where DFT 
predicts a metal). 

Related to the PD is the less known and less worked out use of 2-body wave functions 



2l|. The K 



(geminals) ip K {l,2) for the description of an electronic many-body system 
are 2-body quantum numbers. Recently, in an eventually possible pair density functional 
theory fPDFT') I22I and in the Kimball- Overhauser approach of the homogenous electron 



y,Q,y,HQ 



gas (HEG) 23J, |24|, |25j, |2£|, |27j, |2$, |29j such geminals are discussed as solutions of an effective 
2-body Schrodinger equation and used to parametrize the pair density (PD) p 2 (l,2) = 
Y^k ^K^xiXi 2)| 2 . More general, the so-called natural geminals diagonalize the 2-body 
RDM (2-matrix) 7 2 (1|1', 2|2') = ^^(1, 2)^£(1', 2'). It is quite natural to CI expand 



the ipjc(l,2) in terms of the ip K (l). In Ref. 



30| it is suggested to generalize the Kimball- 



Overhauser approach by using its PD geminals as natural geminals in the spectral resolution 
of 72. This generalization has the advantage that from the natural geminals not only follows 
the PD, but also the 1-matrix 7 using thereby the contraction properties of 72. - Geminals 
appear also in "the antisymmetrized power (AGP) function as a flexible ansatz for fermion 
systems with arbitrary N" and they appear in a generalized DFT in 2-particle space, cf. 
20], p. 325 and refs. therein. Besides, there are links between DMFT and geminal theory 
3Jj. Finally, a general remark of Davidson is "scientists have not yet learned to think in 
terms of (^) geminals rather than in N orbitals" [6j], p. 97. 

In the following, the Kimball-Overhauser approach is summarized from an RDM point 
of view and possible generalizations are presented. 



Basic notation 



The GS energy E of a non-relativistic iV-electron system in Born-Oppenheimer approxi- 
mation is partitioned as E = T[y] +V ex t[p] +Mnt[p2] 3 where the well-known linear functionals 



(using the short-hand notation f dl = f d?ri) 

T[ 7 ] = /<n ( (r lh (l|l')|^, i(ri) = -|-(J-) 2 , 

V«W = /dl (>(l)iw(ri), VW[pa]= /^P2(l,2)i-, e 2 = ^- (1) 

appear. f cx t( r ) is the potential of the nuclei (or the jellium background) binding the elec- 
trons. The 1-matrix 7(1|1'), the electron density p(l) > and the PD p 2 {l,2) > are 
functional of the GS wave function . . .,N) arising from \I>(1, . . . , iV)^*^', . . . , N 1 ) by 
appropriate contractions, which means the operations i' = i and J di. (In the many-body 
perturbation theory this term has another meaning.) 7, p and p 2 follow from the 2- matrix 
(which is defined by the N — 2 contraction of . . . , JV)**(1', . . . , N')), 

7 2 (1|1',2|2') = J ^•_;^ (l,2,3,...)^(l / ,2 , ,3,...), Tr 72 = N(N - 1), (2) 

where the wave-function normalization 

/dl . . . dN 1 < „ . . 

jy| |^(1,---,AT)| 2 = 1 (3) 

is used (each particle configuration is naturally counted only once). 

Contraction sum rules, cumulant expansions, spin structures, and pair densities 

The contraction of the 2-matrix yields 

7(1|1') = ^7 Jd2 l2 (l\l',2\2), p(l) = J D 7 (l|l , )=7(l|l). (4) 
D projects out the off-diagonal elements i! 7^ i. The contraction of p2 = D^ 2 yields 

Tr7 = N follows from Tr72 = iV(iV — 1). p is needed for the expectation value of the 
external (the electrons confining) potential, whereas 7 and p 2 are the quantities to calculate 
the kinetic and the interaction energy, respectively, cf. Eq.((H). 



A fairly natural partitioning of 72 is its cumulant expansion 



72 (l|l',2|2')=7 2 HF (l|l , ,2|2')-x(l|l',2|2'), 
7 2 HF (1|1',2|2') = A 7 (l|l / )7(2|2 / ) = 7(1|1')7(2|2') - 7 (1|2 / )7(2|1 / ). 



(6) 



The index HF means generalized (with non-idempotent occupancies) Hartree-Fock part of 
72. A is the antisymmetrizer. x is the so-called cumulant 2- matrix. (For cumulants, cf. e.g. 
Ref. |3) F° r the diagonal elements, Eq.® is read as 



p 2 (l,2) = p(l)p(2) - 7 (1|2) 7 (2|1) -u(l,2) 
with the cumulant PD u = D%. For \ the contraction sum rule (SR) (JIJ is written as 

I d2 X (l|l',2|2) = ^^(1)^(1 - ^)C(l')- 

re 

This contraction SR contains the normalization SR 

dld2 u{l, 2) = Tr X = Tr 7 (l - 7) = ^ v K {\ - v K 



(7) 



(8) 



Nc. 



(9) 



It defines a quantity c < 1, which is the cumulant PD normalization per particle. This 
quantity is referred to as Lowdin parameter (because Lowdin has asked for the meaning of 
Tr7 2 33]). It vanishes for idempotent 



induced non-idempotency 



14. 
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occupancies and increases with increasing correlation 



exchange u K (1 — u K ), what is called particle-hole symmetry 
The spin structure of 7 is simply 

7(1|1') = o ( W^( r il r 'i)> 7(ri|ri) = pf{r), r = |r 



12 \ ■ Note that Eqs. (jHJ) and (HH) are invariant under the 

3J1. 
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/(o) 



(10) 



where the latter expression holds for a homogeneous system with the dimensionless 1-matrix 
f(r). The property /(0) = 1 makes 7(r 1 |r' 1 ) correctly normalized. Note that for an spin- 
unpolarized system p(l) = p(riV2, besides homogeneity makes p(r) = const = p. 



The spin structure of 72 is 
72(1|1',2|2') = 



35] 



<T1, — <T2 7+l i l i , z \ z ) 



+ o (<Wj<W 2 + $n,o£ ^.oi) 7-(l|l', 2|2') 



(11) 



which gives rise to a singlet 2-matrix 7+ and a triplet 2-matrix 7_ with the normalizations 
Tr7± = |(|±l). For a homogeneous system the diagonal elements of 72 define the 
spin-dependent dimensionless PD g according to p 2 g = p 2 , 

#(1,2) = -(1 - 5 au(72 ) 8 ai _ a2 g + (r 12 ) + ~(1 + 8 ai , a2 ) g_(r 12 ), r 12 = |r x — r a |. (12) 

Here g+(r) is the singlet PD and g~(r) is the triplet PD, which are normalized as 

Jd 3 r |[0 ± ( r )-l] =±1, <7±(oo) = l. (13) 

Thus the PDs for spin-parallel (g p ) and spin-antiparallel (g a ) electron pairs can be defined 
with the help of g + and g- as 



(r) = g-(r), g a {r) = -\g+(r) + g-(r)) 



(14) 



They have the following normalizations 

Jd 3 r |[l-^(r)] = l, Jd 3 r^[l-g a (r)] = 0, g p>a (cx>) = I. (15) 

With these definitions, the spin-summed PD is 

11/* 

9(r) = - [g a (r) + g p (r)) = - [g + {r) + 3g_(r)] ,ld 3 rp[l- g(r)} = 1, g(oo) = 1. (16) 

Note that in Refs.Q, 3 the 9±( r ) are defined differently, namely with a factor |, such that 
we have here <?±(oo) = 1 (whereas in Refs. 2<], 2^ it is \). 

Eqs. ()13|) - (fT6*j) may be equivalently written in terms of the dimensionless cumulant PD h 
defined by p 2 h = u. Thus Eq.(JZJ) takes the form 



g{l,2) = \-\8^\f{r)\ 2 -h{l,2) 



(17) 



from which follows 



(r) = 1 - \f{r)\ 2 - h p (r), g a {r) = 1 - h a (r). 



(18) 



With these definitions the normalization SRs (j!5j) become 



d 3 r ^h p (r) = c, / d 3 r ^h a (r) = 0, h p>a (c^ ) = 0. 



(19) 



Next singlet /triplet terms h± are defined - analog to Eq.(JTH) by h p = h_ and h a = 
Their normalization SRs are 



d 3 r ^rh ± {r) = =pc, 



1-- / d 3 r |/(r)| 2 , /i±(oc) = 0, 

1 



(20) 



what follows from Eq. (jl9j) . Analog with Eq. (|16)l . the spin-summed cumulant PD is 



h(r) = ~[h a {r) + h p {r)} = -[h + (r) + 3/i_(r)], / d 3 rph(r) = c, h(oo) = 0. (21) 
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— ± 1 



(22) 



Similar as for g± the definition of h±(r) differs from the definition used in Refs. 
For the j± of Eq.((TTJ) the contraction SRs 

J d 3 r 2 7 ± (ri|r / 1 ,r 2 |r 2 ) = ^7( r i| r i) h 

hold, what contains the normalization SRs Tr7-|- = y (y ± l) and agrees together with the 
spin structure Eq. fTTJl also with / d27 2 (l|l', 2|2) = j(l\l')(N - 1) of the contraction SR (gj). 
For x±i the cumulant part of j±, defined by j± = 7^ F — x±, it follows from Eq. (j22|) 

1 



/ 



d r2X±(ri|ri,r 2 |r 2 ) = =F; 



7l r il r 



J ^27(^1^)7(^1^) 



(23) 



Trx± = -FyC. This contraction SR does not contain any non-size extensive term contrary 
to the rhs of Eq. (j22j) . The difference 7^ F — j± makes the non-size extensively normalizable 
and contractable terms to cancel each other what is the prerequisite for the thermodynamic 
limit. The contraction SR (|23|) allows one for a homogeneous system to calculate the (non- 
idempotent) momentum distribution n(k) by solving a quadratic equation provided that 
the cumulant matrices x± are known e.g. from perturbation theory [x± is given by linked 
diagrams, but the RPA-like approximation used in Ref. 3] yields only the idempotent 
n°(k) = 0(1 — k)] or from the natural geminals as solutions of an effective 2-body Schrodinger 
equation. What the latter means is described in the following. We start with the 1-matrix 7, 
derive from this 7± F , the HF part of the 2-matrix, and use it in the definition x± — 7± F — 7±- 



Natural orbitals and natural geminals 

For a homogeneous system the natural orbitals of 7(17 Ir^) are plane waves y2k(r) = "^ elkr 
(Q = normalization volume) and their occupancies give the momentum distribution n(k), 



M = I E ^) e " < »(*) < L E = / ( 24 ) 



resulting thus from the Fourier transform of the 1-matrix /(r) 

fid 3 A; 

n(fc) allows one to analyze Compton scattering data. With increasing correlation, the quasi- 
particle weight Zp = ""-(I - ) — n(l + ) (being 1 for 'no interaction' or r s = 0) decreases. 
Note that in the more general spin-polarized case there are two different momentum dis- 
tributions (for spin- up and spin-down). But for 'no polarization' they coincide. For recent 
parametrizations of the momentum distribution(s) cf. 

With this spectral resolution of the 1-matrix 7 one can easily write down 73 = A77 and 
from this follow its singlet /triplet components 7^ F . Thereby the geminals ^ elKR "^ elkr 
with R = 2(1*1 + 1*2) an d r = ri — r2 appear together with the weight n(ki)n(k2) as the 
probability of finding two electron momenta k 12 = |K ± k, where K = k x + k 2 is the 
total momentum and k = |(ki — k 2 ) is the half relative momentum. The first factor 
describes the free-particle center-of-mass motion, whereas ^e lkr = 4= i'j£,(fcr)Y£ (e^), 
ji{kr) = ji(kr)Yi(e r ), L = (I, mi) describes the relative motion for 7^. Starting with 
7^ F = A77 and defining 7^ F according to the spin structure Eq. fjllj) . it results 

7 f = ^^i)«(^)v±( r i> r 2;ki,k 2 ) ip^^-MM) (25) 

ki,2 

with ipl = ^e iKR ^ i l j L (kv)Yl(e k ) and 

Tr 7 « F = Y,<h)n{k 2 ) [1 ± 5 kllk2 ] = y (j ± l) T f c 



(26) 



So, the spectral resolution (|23|) can be written in terms of free-electron geminals ji(kr) and 
an occupancy matrix 

/TO 
-^■Yl{e k )n{k 1 )n{k 2 )Y I/ (e k ). (27) 

Its diagonalization yields an L-mixing j^. 

With the aim to get finally x± = 7± F — 7±> h is assumed that 7± has the same form as 
7^ F in Eq. (|25jl with the only difference of replacing the free-electron geminals (p°± or ji(kr) 



by interacting-electron geminals (p± or Ri(r, k): 

7±(rir' 1 , r 2 |r' 2 ) = ^ n(k 1 )n(k 2 )ifi±(r 1 , r 2 ; k x , k 2 ) <£>±( r 'n r 2, k i, k 2)- (28) 

ki, 2 

If these interacting-electron geminals Ri together with the momentum distribution n(k) are 
available, then the PDs g± follow from 

± 2 

— V „,(k) R?(r. k). u(k\ = 

N 



L k 



(29) 



K 



with fEkM^) = 1 an d MO) = 2 3 (1 



Assuming that the Ri are scattering solutions of a radial Schrodinger equation (cf. next 
Sec.) with a large-r asymptotics (phase shifted compared with j t ) according to 

1 



1 7T 

Ri(r, k) -> — sin (At - /- + ^(A;)), jj(Ar) 



the normalization SRs 1)1 Hj) or ()2()jl take the form 



.7T, 



, ■ sin(Ar — /— ), 
kr v 2 7 ' 



(30) 



7T 



E 



dk 



r)i(k) = ±c. 



(31) 



One may compare these SRs with the well-known Friedel SR for point defects in metals with 
their screening cloud around an impurity. If the spectral resolutions ()25|) and ()28|) are used 
in x± — 7± F — 7±, the contraction SRs (|23|) can be written as [^J 
^ ± 



7T 



E 



poo 

/ rfA; 


d2 3 n(2K,k)~ 


dk 


/o 



7]{k) + 6±(«) = ±7l(«)[l - 7l(K)], 



//(K,A;) 



4tt 



■n(k 1 )n(k 2 ). 



(32) 



The quantity 6±(k) is defined in the next Section. With X^k,^±( k ) = 0> the contraction 
SRs (j3~2*j) contain the normalization SRs (|3*T|) as special cases. Thereby 4 2 3 /z(2k, fc) = 
/i(fc) is used. The contraction SR (|3*2)l may be considered as the spectral resolution of the 
normalization SR ()31j1 . Both Eq. ()31j1 and Eq. ()32)1 are relations between the 1-body quantity 
n(k) and the 2-body quantities Ri(r, k) in addition to the virial theorem. 



The Kimball- Overhauser approach and its generalizations 

Now the question is, where to get from the geminals Ri(r,k). In Refs. 0, ^ 
is intuitevely assumed that they are (at least for the parametrization (|29|) of the PDs) the 
solution of the radial Schrodinger equation 

Ri(r,k) = 0, (33) 



1 d 2 1(1 + 1) 
—£-»t H « \-v±{r) 



k 



which arises from a 2-body Schrodinger equation with an effective interaction potential 



v±(r) 



v scr( r ) possibly different for + (= even I) and - (= odd /). This is the Coulomb 



repulsion between two electrons effectively screened by the Fermi-Coulomb hole around each 
electron. This therefore attractive screening potential comes in the Hartree description 2(| 
from the Poisson equation Av scr (r) = 47rp[l — g(r)], which makes the approach a self- 
consistent one: g°(r) — > v±(r) — > R®(r,k) — > g 1 ^) —>■■■■. The results of this approach 
are in excellent agreement with the quantum-Monte-Carlo data of Refs. fl 0- They 
have been further improved with the inclusion of exchange and correlation in Ref. 



Meanwhile also the spin-polarized HEG has been treated in this way |2a]. These successes 
say that this approach contains at least some truth, even if it is not exact and it confirms 
the above assumption that only the wave functions change from (f° to if when going from 
the known 7± F to the unknown 7± leaving the occupancy weight unchanged. 

One generalization is the assumption that the PD geminals of Eq. (}33|) can be used in the 
spectral resolution (}2*Hj) also as 2-matrix geminals. This allows one to calculate not only the 
PDs g±(r) but also the momentum distribution n(k) through the contraction SRs ()32)1 - at 
least in principle. The quantity b±(r) therein is defined by 



1 ± 1 



N 2 

L,L' K12 



X 



2 x - 

[hlv(2(k - Kt), k) - im w (2(k - K 2 ),k)] 



(34) 



dk 2 dk 2 

where v±(ku) and -Rl(k, k) are the Fourier transforms of v±(r) and Rl(y, k), respectively. If 
one starts with a PD g°(r), then it follows (e.g. in the Hartree approximation) the effective 



interaction potential v±(r) yielding geminals R^(r,k) with their phase shifts t]f(k) and a 
PD g x {r). Next with a starting momentum distribution n°(k) the lhs of Eq. ((321) can be 
calculated. The result is a quadratic equation. Its solution gives n 1 ^), etc., until self- 
consistency is reached finally. Whether this really works has to be checked. For the "PD 
to 2-matrix" generalization in terms of phase shift SRs for normalization and contraction 
cf. Refs. Q and H, 

respectively. For the generalization of these SRs to the case of the 
spin-polarized HEG cf. Ref. 35^ . 

Another generalization concerns inhomogeneous systems. It has been already discussed 



in Ref. 



26[. Here the problem is alternatively viewed for the case of an extended system 



(e.g. a jellium with a crystalline periodic background density). Then one has to solve the 
2-body Schrodinger equation 



I [t{ri) + iwfo) + v n (n)] + v±{r h r 2 ) - + k 2 2 ) 

U=l,2 



<p ± = (35) 



with the local Hartree potential v#(r) and an appropriately screened Coulomb repulsion 
v±(ti,T2). The self-consistent procedure would have to start with a reasonable approxima- 
tion for the natural orbitals and their occupancies. ^From this follows 7± F . From the solution 
of Eq. (jHHjl [where on the Hartree level, £2(1*1, r 2 ) is needed as an input for v±(tx, r 2 )], follows 
also 7-j-. This gives the new PD. But also x± = 7± F — 7± is then available, which - used in 
Eq. (jEJ - yields new natural orbitals and occupancies. 



Summary and outlook 

The Kimball-Overhauser approach for the pair density of the spin-unpolarized homoge- 
neous electron gas (HEG) in terms of geminals is revisited from a reduced-density-matrix 
(RDM) point of view. We start with the definition of the 2-body RDM and its spectral 
resolution in terms of natural geminals. This gives quite naturally the Kimball-Overhauser 
parametrization of the pair density in terms of such geminals. Thereby it is assumed that 
the geminals which parametrize the pair density can be used also as natural geminals. This 
assumption has the advantage that also the 1-matrix can be calculated as a consequence of 



the contraction propeties of the 2-matrix. An important role plays the size-extensivity of 
the cumulants, which allows one to consider the thermodynamic limit. Open questions are: 

• Does the asymptotics Eq. lfHU]) of Ri(r,k) give [via Eq. (J2HJ)] the correct asymptotics 
of g±(r)7 If this is not the case, can this deficiency be removed by using non-local 
effective interaction potentials ? 

• Why the lhs of Eqs. (jHljl and (|H2jl are note particle-hole symmetric? 

• Is there a link to the concept of strongly orthogonal geminals? 

• What are the contraction sum rules for the spin-polarized HEG j^|? 

• Can the effective 2-body Schrodinger equations (}33|) . (J35|) be derived from the hierarchy 
of contracted Schrodinger equations or from the Bethe-Salpeter equation or • • • ? To 
what extend is the Kimball- Overhauser approach and its generalizations related to 



a possible pair-density functional theory 22 1 and to 2-body cluster expansions 39(? 
How to treat finite systems in terms of a screened Coulomb repulsion? 

What are the peculiarities of x± an d v± for bondbreaking situations, for metals, semi- 
conductors, ferromagnets, antiferromagnets, superconductors, off-diagonal-long-range 
order, ferromagnetic superconductors, mixed valence compounds, heavy-fermion sys- 
tems, non-Fermi liquid behavior, quantum criticality, • • ■ ? 



In Ref. [40] calculational electronic-structure methods are reviewed including attempts 
to generalize density-functional theories (e.g. density-matrix functional theory), to further 
develop the RDM theory (contracted Schrodinger equations, iV-represent ability to 
make the accurate quantum-chemical methods (configurational interaction, coupled cluster, 
M0ller-Plesset) applicable also to extended systems (crystalline solids), and to treat systems 
with strong electron correlations from first principles. 
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